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SUMMARY 

The steady state displacements of a rotating advanced turboprop, SR3 , were 
computed using the geometrically nonlinear capabilities of COSMIC NASTRAN 
Rigid Format 4 and MSC NASTRAN Solution 64. Displacements were computed for 
rotational speeds up to 10 000 rpm. The results showed that a complete non- 
linear analysis Is required. Neither a linear elastic analysis nor a one step 
differential stiffness analysis were adequate. The Inaccuracy of these linear 
analyses Increased with Increases In rotational speeds. 

A modified Newton-Raphson algorithm used by MSC NASTRAN Solution 64 and an 
Iterative scheme used by COSMIC NASTRAN In Rigid Format 4 were employed for 
predicting geometrically nonlinear response of SR3 . Solution 64 was run two 
ways: (1) using a constant centrifugal body force computed from the undeform- 

ed blade, and (2) using a force field corrected In each Iteration for changes 
In the blade's position. Rigid Format 4 does not possess this load updating 
capability. When the load updating feature of Solution 64 was not used, the 
two programs generated similar tip deflections and rotations. When the load 
updating capability of Solution 64 was used, the tip displacement varied sig- 
nificantly between MSC Solution 64 and COSMIC Rigid Format 4. 

The SR3 turboprop blade has been shown to behave In a geometrically non- 
linear fashion under centrifugal loading. This nonlinear response becomes 
Increasingly Important at higher rotational speeds. Also, since MSC NASTRAN 
Solution 64 can correct the centrifugal load for changes In blade position. It 
Is a more desirable analysis tool for predicting nonlinear behavior than 
COSMIC NASTRAN Rigid Format 4. 


INTRODUCTION 

Increasing concern with Improving aircraft fuel efficiency has brought 
about renewed Interest In propeller propulsion systems. Improved multl-bladed 
propellers, termed advanced turboprops, have the potential for reduced fuel 
consumption while maintaining the performance levels of modern turbofans. 

These turboprops feature thin, flexible, swept blades with complex structural 
properties. Extensive research In both analytical and experimental techniques 
has been and continues to be conducted by researchers In the field to better 
understand the structural and aerodynamic response of these complicated blades. 

An Important area to Investigate Is the computation of steady state de- 
flections In rotating turboprop blades. Previous experience has demonstrated 
that linear analysis techniques are Inaccurate for predicting the deflections 



of these rotating blades. Because the blades are relatively flexible they 
respond to loading with relatively large deflections. This complicates the 
calculation of the blade's stiffness since the stiffness and deflections are 
mutually dependent. The computation of this behavior requires a geometrically 
nonlinear analysis. 

The analysis Is further complicated by the fact that the centrifugal loads 
are also displacement dependent. Since centrifugal force Is proportional to 
the radius from the rotational axis the magnitude and distribution of cen- 
trifugal loads will change as the blade displaces. The research presented In 
this paper Investigates this effect. 

The purpose of the research presented In this paper Is to Investigate the 
use of NASTRAN, a well-known finite element program, for predicting steady 
state deflections of advanced turboprops subject to centrifugal loading. This 
computer program was selected due to Its 1 nonlinear analysis capability and 
availability. Furthermore, It Is universally accepted and has been used suc- 
cessfully to analyze geometrically nonlinear structural problems In the past. 
Two versions of NASTRAN were employed, the first was MSC (Macneal- 
Schwendler Corporation) NASTRAN Solution 64; the second was COSMIC (Computer 
Software Management and Information Center) NASTRAN Rigid Format 4. Both of 
these programs can be used for geometrically nonlinear analyses. It Is con- 
venient to use COSMIC NASTRAN Rigid Format 4 because the stiffness matrix 
generated In this program Is In a format compatible for Input Into the sub- 
sequent aerodynamic flutter analysis. The turboprop blade was analyzed using 
MSC NASTRAN Solution 64 In order that the results from rigid format 4 could be 
compared to an Independent nonlinear analysis. 

A representative advanced turboprop, named SR3, was used as the model for 
Investigating the use of NASTRAN for computing deflections of advanced turbo- 
props. SR3 has most of the structurally related characteristics found In 
typical advanced turboprop blades (fig. 1). 


BACKGROUND 

At NASA Lewis Research Center COSMIC NASTRAN runs on a UNIVAC 1100 com- 
puter system which utilizes a 32 bit processor and uses double precision. MSC 
NASTRAN Is run on a Cray 1-S. This program uses single precision and a 64 bit 
processor. 

As previously mentioned. It Is desirable to use COSMIC NASTRAN for the 
nonlinear analysis because the results from this program can easily be used In 
the aerodynamic flutter analysis. Since both the nonlinear analysis (Rigid 
Format 4) and the aerodynamic flutter analysis (Rigid Format 9) are run on the 
UNIVAC using COSMIC NASTRAN, It Is relatively simple to transfer results from 
one program to the other. More specifically. It Is necessary that the stif- 
fness matrix generated from the nonlinear analysis In Rigid Format 4 be trans- 
ferred to Rigid Format 9 for use In the flutter analysis. If the flutter 
analysis does not use a stiffness matrix that Includes nonlinear effects, cor- 
rect frequencies will not be computed and the results of the flutter analysis 
will be In error. The stiffness matrix from MSC Solution 64 cannot be used 
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for the flutter analysis because there Is presently no practical method for 
transferring matrices between MSC and COSMIC NASTRAN. An effort Is currently 
being made to eliminate this problem. 

From among all the MSC and COSMIC NASTRAN solution sequences available for 
geometrically nonlinear response Solution 64 Is thought to be the most 
complete. In order to verify the accuracy of Rigid Format 4, the results from 
this program are compared to the results from MSC NASTRAN Solution 64. 


MSC NASTRAN Solution 64 for Geometrically Nonlinear Analysis 

MSC NASTRAN Solution 64 provides a straightforward means for performing 
geometrically nonlinear analysis utilizing a modified Newton-Raphson algorithm 
(refs. 1 and 2). The algorithm used In this solution sequence Is described In 
the MSC NASTRAN application manual and will be repeated here for the reader's 
convenience. To best understand the Newton-Raphson algorithm used by MSC 
NASTRAN, a simplified Incremental scheme will first be described and then It 
will be shown how this scheme may be expanded Into the Newton-Raphson 
algorithm. 

The objective of a nonlinear analysis Is to simulate the correct displace- 
ment versus load relationship. One means of accomplishing this Is to use an 
Incremental algorithm In which the nonlinear response Is divided Into a series 
of linear steps, each step representing an Increase In load (fig. 2). 
Displacements are accumulated after each step and a new tangent stiffness Is 
calculated based on the structure's deformed position. By using this approach 
the nonlinear response can be approximated by a series of linear segments. 

The load Increment determines how close the approximation will match the 
actual response curve. The larger the load Increment, the further the computed 
response will be from the actual response. When large load steps are used, 
the purely Incremental algorithm will deviate considerably from the correct 
structural behavior. As shown In figure 2, while the actual displacement for 
applied load R Is at D a , the computed displacement Is only at Do. 

This deviation results In equilibrium conditions not being satisfied at the 
end of the load step. While the difference AR between externally applied 
loads and Internal element forces should equal zero, they are actually equal 
to a nonzero value. 

The error produced by the Incremental approach can be eliminated In the 
limit by requiring that equilibrium conditions be satisfied at the end of each 
load step. This approach Is analogous to the conventional Newton-Raphson 
algorithm. The MSC NASTRAN form of the Newton-Raphson approach requires that 
equilibrium be satisfied (AR = 0) by balancing the externally applied loads 
and the Internal forces at the end of each load step. This constraint pro- 
hibits the computed response from deviating from the actual response. This 
algorithm has the added advantage that, for typical structures where the load- 
displacement curve Is relatively smooth, the entire load can be applied In the 
first step without effecting the accuracy of the final results. 

The MSC NASTRAN form of the Newton-Raphson algorithm Is based on an Iter- 
ative solution of the equation of equilibrium. In matrix notation, this 
equation takes the form: 
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[Ki] {AD i+1 } = {ARiJ 


( 1 ) 


-x/- x 


Where the load Imbalance {aR]} = { R a } - Z[k]] {d ^ } . [K] ] and [ k ^ ] 
are the global and element stiffnesses respectively, {R a } Is the applied 
load, {d]} are the element displacements, and {AD] + ]} Is the Increment 
In the global displacement at the end of the Iteration. The global stiffness 
matrix [K]] Includes both the elastic and differential stiffness matrices. 


The Newton-Raphson algorithm Is shown graphically In figure 3. For the 
Initial Iteration, the structure begins In Its undeformed position with the 
Internal forces set to zero (z[k]] {d]} = {0}. This reduces the right side 
of equation (1) to the applied load vector {R a }. 

By using the load {R a } and the tangent stiffness [K 0 ] at the origin, the 
displacement {AD]} Is computed. After applying {A0-|} to the undeformed 
blade, the position of the displaced blade Is established. Using this deformed 
shape, the Internal forces (z[k]] {d^}) are generated along with a new 
tangent stiffness. The next Iteration Is Initiated by computing a new set 
of displacements using the Imbalance ( { R a } - z[k]] {d]}) as the applied 
load and the new tangent stiffness as the stiffness. This procedure Is re- 
peated until the load Imbalance Is reduced to within the desired tolerance. 
Reducing the tolerance produces more accurate results but requires additional 
Iterations . 


It Is not mandatory for the load vector { R a } to be equal to the 
total external load In every Iteration. Instead, { R a } can be only a 
fraction of the total external load In early Iterations and then can be In- 
cremented to the full load In the final Iterations. In the example shown In 
figure 3, {R a } Is taken as the full externally applied load In every 
Iteration. If the structure Is highly nonlinear and the total external load 
Is applied In the Initial Iterations the structure may move to an unstable 
condition. The advantage of using load Increments Is that the path to 
Instability can be monitored. 

In the modified Newton-Raphson approach the elastic stiffness matrix Is 
not updated In each Iteration. Instead, the original stiffness matrix derived 
from the undeformed geometry Is used for every Iteration. Using this original 
stiffness does not affect the final results because the final displacement Is 
Independent of the path used! (The element stiffnesses which are used to 
compute the load Imbalance on the right side of eq. (1) are updated In each 
Iteration.) By referring to figure 3, one can visualize the affect of using 
the same stiffness In every Iteration. If [Kq] Is used as the stiffness in 
every Iteration, the Incremental displacements will be different from the 
displacements shown In figure 3 (altering the path), but the final converged 
displacement will still be equal to D a . 

The advantage of using an unaltered stiffness matrix Is that the cost of 
forming and decomposing a new elastic stiffness matrix for each cycle Is 
eliminated. The disadvantage Is that the solution will require more iter- 
ations to converge. Although the elastic stiffness need not be updated with 
MSC NASTRAN, the user does have the option of computing a new differential 
stiffness matrix In any Iteration. It Is useful to be able to update the dif- 
ferential stiffnesses matrix since It Is a function of element stresses and 
can change considerably with load Increments. 
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MSC NASTRAN Solution 64 uses "subcases" to control the execution of the 
nonlinear analysis. Specifically, the number of subcases specified In the 
NASTRAN Input data deck determines the number of Iterations that will be per- 
formed. The user must specify at least two subcases. In the first subcase a 
linear elastic analysis Is used to compute an Initial deflected shape. This 
displaced shape Is then used In the second subcase to compute the differential 
stiffness matrix along with a new set of displacements. Subsequent subcases 
are used for Iterating on the equilibrium equation (eq. 1). 

Solution 64 has two Important features that Increase the flexibility of 
the nonlinear analysis. First, Solution 64 has the, ability to recompute the 
external loads before each Iteration. This Is particularly advantageous for 
centrifugal loads since they are dependent on the radius of the blade's mass 
from the rotational axis and will, therefore, change In each Iteration as the 
blade deflects. The second asset of Solution 64 Is Its ability to apply a 
fraction of the full external load In earlier subcases. This advantage has 
previously been discussed. Neither of these capabilities are available In 
COSMIC NASTRAN Rigid Format 4. 


COSMIC NASTRAN Rigid Format 4 for Differential Stiffness 

Rigid Format 4 of COSMIC NASTRAN Is designed to solve geometrically non- 
linear problems (ref. 3). This rigid format uses an Iterative solution 
sequence based on an extension of the one step differential algorithm. The 
governing equation for this rigid format Is: 

[[K] + [K d (Di)]J {D U1 } = { R a } (2) 

Where [K] Is the global elastic stiffness, [K d ( 0^ ) ] Is the differential 
stiffness computed for the displaced blade at { D^ > , { R a > Is the 
applied load, and Is the displacement computed at the end of the 

Iteration. In this equation the differential stiffness matrix [K d ] Is up- 
dated In each Iteration whereas the elastic stiffness [K] remains constant. 
Equation (2) Is not guaranteed to converge to the correct displacement because 
this equation, unlike equation (1), contains no check to Insure a small load 
Imbalance between Internal element forces and external body forces. Thus, for 
highly nonlinear structures, equation (2) may not converge. 

Rigid Format 4 requires two subcases. The first subcase Is used for per- 
forming linear elastic analyses without Incorporating differential stiffness. 
The resulting element stresses from the first subcase are used In the second 
subcase to compute the differential stiffness matrix. Equation (2) Is then 
repeated until the weighted difference between the differential stiffness In 
subsequent Iterations Is within the desired tolerance or until the maximum 
desired number of Iterations has been completed. Both the tolerance and 
number of Iterations can be controlled by the user. The linear elastic and 
Initial differential analyses performed by Rigid Format 4 are comparable to 
the first two subcases of MSC NASTRAN Solution 64. 

The major limitation of using Rigid Format 4 for computing deflections of 
rotating structures Is the Inability to compensate for the change In the 
centrifugal body force as the blade displaces. This limitation can Introduce 
considerable error. 
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RESULTS 


A comparison of results between MSC NASTRAN Solution 64 and COSMIC NASTRAN 
Rigid Format 4 was made using the advanced turboprop model SR-3. This turbo- 
prop blade Is a small scale titanium model measuring 12-1/4 In. from the 
rotational axis to the blade tip. The blade thickness decreases from over one 
In. at the root to 0.016 In. at the tip. As shown In figure 1, the blade Is 
both highly swept and twisted. This geometry Is Intended to optimize aero- 
dynamic performance and minimize noise while maintaining structural Integri- 
ty. The finite element model consists of 346 triangular plate elements and 
206 grid points. The "CTRIA2" and "CTRIA3" elements were used for the COSMIC 
and MSC runs, respectively. The plate element formulations Include both 
membrane and bending action. The base of the root of the blade Is modeled as 
fully constrained. 

MSC NASTRAN Solution 64 was used to compute steady state displacements of 
SR-3 at 3500, 7000, 8600, and 10 000 rpm. Two sets of displacement data were 
generated. The first set was generated using a constant centrifugal body 
force computed from the original, undeformed blade. The second set used a 
body force updated In each Iteration based on the position of the deformed 
blade. Seven subcases (5 Iterations) were used for Solution 64. Convergence 
was evaluated by comparing the difference In displacements between the sixth 
and seventh subcases which were on the order of only one tenth of one percent. 
The five Iterations produced a ratio of unbalanced forces at unconstrained 
nodes to forces at the constrained root of 10 -4 . This small ratio Indicates 
that the load Imbalance at the unconstrained nodes Is relatively small and 
that equation (1) has adequately converged. 

Steady state displacements were computed using COSMIC NASTRAN Rigid Format 
4 at the same rotational speeds as MSC Solution 64. Only two Iterations are 
required to produce converged displacements at 3500 rpm. Additional Iter- 
ations changed the final displacements by only 1/2 percent. The number of 
Iterations was Increased to three by decreasing the default tolerance from 
10“5 to 10“6 (PARAM, EPSIO, l.E-6). At 10 000 rpm seven Iterations were 
required for the same order of accuracy. 

Similar finite element formulations are used for both COSMIC NASTRAN and 
MSC NASTRAN. As expected, the two programs computed similar, though not 
identical, displacements In the linear elastic subcase. The small discrepancy 
In the displacements Is probably the result of either the minor difference In 
element formulations or the difference In criteria used by the two programs 
for testing small values of stiffness In the global stiffness matrix. Since 
the formulation for triangular plate elements does not Include In-plane 
rotational stiffness, the assembled global stiffness matrix may contain very 
small or zero entries. These small values are usually removed by constraining 
appropriate rotational degrees of freedom. The default minimum value used by 
COSMIC NASTRAN Is larger than the value used by MSC NASTRAN. Thus more con- 
straints are present In the COSMIC runs. These additional constraints add 
stiffness to the blade and typically lower the Rigid Format 4 elastic dis- 
placements. However, the displacements In the subsequent nonlinear analysis 
do not appear to be affected by the additional constraints. 

Tip deflection (defined as total deflection of midchord at blade tip) as a 
function of rotational speed Is shown In figure 4. The general trend Is for 
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deflections to Increase up to around 7000 rpm and then level off. All three 
analyses follow this trend. COSMIC Rigid Format 4 and MSC Solution 64 
produced similar results when centrifugal loads were not updated with changes 
In blade position. The difference In tip deflections between these two 
solutions with no load update decreased from 14 percent at 3500 rpm to only 5 
percent at 10 000 rpm. There Is a constant 0.01 Inch difference In displace- 
ment between the two solutions above 3500 rpm. When the load updating 
capability of Solution 64 was utilized, the deflections were considerably d 1 f 
ferent. The deflection curve generated using the load updating algorithm 
shows larger deflections than the other two curves (32 percent larger at 
10 000 rpm). This was expected because the centrifugal force actually In- 
creases as the blade straightens out and the radius of the blade's mass from 
the rotational axis Increases. It Is evident that at high rotational speeds a 
load updating algorithm must be used to account for the large changes In blade 
position. 

Figure 5 shows blade tip rotation as a function of rotational speed. Tip 
rotation Is defined as blade tip chord twist about the pitch axis In a plane 
normal to the pitch axis. The pitch axis Is shown In figure 1. The general 
trend In this figure Is for the blade to untwist as rotational speed Increas- 
es. As with tip deflections, Solution 64 with load updating produced the 
largest rotations. The differences between the load updating and constant 
load solutions were not as notable for tip rotations as they were for tip mid 
chprd deflections. The variation In tip rotations between Rigid Format 4 and 
Solution 64 with load updating was 19 percent at 10 000 rpm. 

Figure 6 shows a comparison between the three levels of analysis used In 
MSC NASTRAN Solution 64. As expected, the linear elastic analysis (subcase 1) 
over predicted the tip deflection. This Is due to the stiffening effect from 
the centrifugal force field which was not Included In the linear analysis. 

The relationship between deflection and rotational speed Is parabolic since 
centrifugal force Is proportional to the square of rotational speed. The 
deflection computed In the linear analysis Is In contradiction to the trend 
predicted by the nonlinear analysis shown In the same figure. The nonlinear 
analysis shows the tip deflection leveling off. The displacements from sub- 
case 1 are considerably greater than those computed In subcase 2, which 
Includes centrifugal stiffening effects. While the tip displacement was over 
0.80 In. at 10 000 for the linear solution, It was only 0.21 In. for the line- 
ar elastic plus differential analysis. The difference In deflections between 
the linear and nonlinear solutions Is shown by comparing subcases 1 and 2 to 
subcase 7. The difference In tip deflection at 10 000 rpm, between subcase 7 
and subcase 1 and 2 Is 75 and 33 percent, respectively. This large variation 
between the nonlinear and linear subcases Is a good Indicator of the blade's 
strong nonlinear behavior. 


CONCLUDING REMARKS 

Because of Its ability to update the displacement dependent centrifugal 
force during the solution process, MSC NASTRAN Solution 64 was shown to be 
superior to COSMIC NASTRAN Rigid Format 4 for the geometrically nonlinear 
analysis of advanced propeller blades. However, It cannot be concluded that 
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MSC NASTRAN Solution 64 accurately predicts the steady displacement of advanced 
turboprop blades until additional analytical and experimental studies are com- 
pleted. 
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